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An evolutionary model for the gamma-ray system 
PSR J1311-3430 and its companion 
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ABSTRACT 

The most recent member of the milhsecond p ulsar with very low- m ass companions 
and short orbital periods class, PSR J1311-3430 (jPletsch et al.ll2013 ) is a remarkable 
object in various senses. Besides being the first discovered in gamma-rays, its measured 
features include the very low or absent hydrogen content. We show in this Letter that 
this important piece of information leads to a very restricted range of initial periods 
for a given donor mass. For that purpose, we calculate in detail the evolution of the 
binary system self-consistently, including mass transfer and evaporation, finding the 
features of the new evolutionary path leading to the observed configuration. It is also 
important to remark that the detailed evolutionary history of the system naturally 
leads to a high final pulsar mass, as it seems to be demanded by observations. 

Key words: (stars:) pulsars: general — (stars:) pulsars: individual (PSR J1311-3430) 
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1 INTRODUCTION 

, The discovery of the 2ms, 93 min orbit, "black widow" (BW ) 
■ pulsar PSR J1311-3430 in gamma-rays (|Pletsch et al.ll2012l ) 
' comes to increase the number of such systems in a very inter- 
esting fashion. This is the first time that gamma-ray pulsa- 
tions led to an identification withou t the need of low-en ergy 
counterparts, and additional work (|Romani et al.ll2012h in- 
dicated the absence or very low level of H in the companion, 
plus the suggestive evidence of a large mass of the pulsar 
itself. These facts are quite important to reconstruct the 

evolutionary path that led to the present state. 

In a former work l|Benvenuto. De Vito. fc HorvathI 

[2013) we have modelled a similar system, namely the bi- 
nary millisecond pulsar PSR J1719-14 38 in a 2.2 h orbi t 
featuring a very low mass companion (Bailes et al. 201lJ). 
We shall show now that the additional information gathered 
from PSR J1311-3430 is useful to refine the models and re- 
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strict the possible initial states to a rather small range, for 
a given assumed value of the initial donor mass. 

We have described the features and tests of the 

evolutiona ry code employed in the se studies elsewhere 

(see Beiivenuto. De Vito. fc HorvathI |2012| : and references 
therein). A simultaneous integration of the stellar structure, 
mass transfer rate (if any) and orbit evolution is performed 
taking into account i) accretion onto the neutron star (NS); 
ii) the evaporating wind driven by the pulsar radiation; and 
iii) the irradiation feedback, present when the donor star 
that transfers mass onto the NS and accr etion illuminates 
the d onor star modifying its evolution fsee lBiining fc Ritterl 
I2OO4I for further details on this item). 

The three effects have been shown to be important for 
the evolution of PSR J1719-1438. The case of PSR J1311- 
3430 adds considerable new information to address because 
spectro scopic/photometric data has shown that H is almost 
absent (jRomani et al.ll2012l ). at the level of a number abun- 
dance uh < 10~^ in contrast with other members of the 
groufij; and also that a high-mass for the pulsar is favoured, 
depending somewhat on the interpretation of the light curve 
but bounded by MpsR ^ 2.1 Mq. We shall show now that 



^ It corresponds to a H mass abundance X < 2.5 X 10 ® for a 
helium-dominated composition. 
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the evolutionary models can match these and the rest of the 
observed parameters, constraining quite tightly the initial 
state of the system. 



2 CALCULATIONS 

The evolution of PSR J1311-3430 system has been 
modelled within the same scen ario depicted in 
iBenvenuto. De Vito. fc HorvathI l|2012r ) starting with a 
normal, solar composition star and a just-formed NS 
orbiting together with an initial period P; < 1 d. 

When the donor star radius R2 equals the radius of 
the Roche Lobe -RzO, it occurs the onset of the Roche Lobe 
Overflow (RLOF). Around this epoch tidal dissipation forces 
quickly circularise the orbit. The donor star transfers mass 
across the Lagrangian point Li towards the NS, causing the 
orbit to evolve. As in IPodsiadlowski. Rappaport. fc Pfahll 
l|2002l) . we have parametrized the uncertain fraction of mass 
effectively accreted by the NS with a quantity /3 > (that 
is, Afi — —J3M2), assuming that it is always below the Ed- 
dington limit Msdd — 2 x 10~* Mq yr~^. Since, in general, 
/? < 1 some material is lost from the system, carrying away 
specific angular momentum of the secondary. Fortunately, 
the value of /3 is not c ritical in determining the e volution 
of this kind of systems l|De Vito fc Benvenutdl20ia ). and so 
we have as sumed an average value o f /3 = 1/2. Gravitational 
radiation (ILandau fc Lifshita [l97l|) and magnetic braking 
l|Verbunt fc Z waan 1981) are known to provide relevant an- 
gular momentum sinks. Additionally, the M2 due to RLO F 
has been described by the e xpressions given in Rittei 1119881). 

As stated in IBenvenuto. De Vito. fc HorvathI 

20121). this "standard" prescripti on valid for LMXBs 

Podsiadlowski. Rappaport. fc Pfahll I2OO2I ) must be 

supplemented with a wind evaporation Mi.evap law. 

Here, we shall assume the exp ressions given by 

IStevens. Rees. fc Podsiadlowskil (Il992l ). namely 



M: 



2,evap — 






(1) 



where the spin down luminosity Lp of the pulsar is given 
by Lp — A-K^ IiPiPi (Ji is the moment of inertia of the 
NS, Pi is its period and Pi its spin-down rate), U2,esc is 
the escape velocity from the donor star surface, a is the 
semi- axis of the circular orbit and / is an effi ciency factor. 
As in IBenvenuto. De Vito. fc HorvathI (|2012r ). we shall set 
/ Lp = 0.04Lq. While irradiation feedback is a fundamental 
ingredient in computing the mass transfer evolution at early 
stages, we have found that it does not affect the evolutionary 
path of the system in the mass-orbital period plane. As we 
are interested here on the final evolution of the system, in 
this Letter we shall neglect irradiation feedback. 

The evolution of the system has been computed 
wi th our detailed (Henyey) e voluti onary code described 
in IBenvenuto fc De Vit"ol l|2003l ) and iDe Vito fc Benvenuto! 
l|2012r ). We started with a fiducial donor mass of 2 Mq to- 
gether with a 1.4 Mq NS, as in the case of PSR J1719-1438. 
The range of initial donor masses is bounded from below 



^ We shall refer the NS (donor star) as the primary (secondary), 
with subindex 1 (2). 



by the isolated-star evolution timescale (M2 ~ 1.0 Mq), 
which should be short enough to allow for the onset of 
mass transfer; and from above by M2 ^ 3.5 Mq, be- 
cause above that initial value the mass tra nsfer is unstable 
jPodsiadlowski. Rappaport. fc Pfahll I2OO2I ). As in the case 
of the evolution of PSR J1719-1438, P must be very short. 
However, extremely short Pi would cause the mass transfer 
to start at the ZAMS; while too long Pi values would ren- 
der a detached wide orbit after a few Gyr, not the observed 
BW-type system. 

Even if these conditions are essentially the same as the 
ones we studied for PSR J1 719-1438, the strong upper limit 
set bv iRomani et al.l (|2012l ) on the H abundance allows an 
important refinement of the models. Its solution rests on a 
delicate interplay between the dynamical evolution of the 
close binary system (CBS) and the structure of the donor 
star in a novel fashion, described as follows. 

If the system started at low periods, around Pi = 0.7 d, 
mass transfer starts when H core abundance is X^^ ~ 
0.39. From this initial condition on, Tc drops slowing down 
H-burning. When M2 « 0.053 Mq the star becomes com- 
pletely convective making it chemically homogeneous up to 
the photosphere with X = 0.099. Since for a larger Pi the 
onset of the RLOF occurs later, X^^ will have a lower 
value and X gets smaller. Eventually, there is a Pi value 
P « 0.86 d for which X^'^"^ = 0, this condition holding 
for every period longer than that. 

For this kind of CBSs, there exists a bifurcation period 
value Pb'. if Pi > Pt the CBS evolves to an open configuration 
forming a low-mass helium white dwarf, while Pi < Pj, leads 
to the kind of BW systems we are interested in. For Pi < 
Pb, the minimum orbital period attained during evolution 
decreases when Pi increaseqj. This behaviour can be seen in 
Fig. [1] where for the system considered here the bifurcation 
period is between 0.880 and 0.890 d. 

For Pi < Pb there occur conditions for which X^^ — 
and mass transfer is able to remove most of H-rich outer 
layers. The donor star becomes a very low-mass He star. 
Remarkably, such configuration leads to masses and orbital 
period in excellent agreement with observations. This tight 
interval 0.87 ^ P/d ^ 0.90 for PSR J1311-3430 is the only 
interval accepted in these calculations, a shorter Pi would 
produce H-rich BW systems. These features are clearly seen 
in Fig. d 

We woul d like to emphasise that a large value of the 
pulsar mass (|Romani et al.ll2012l ) is a natural outcome of 
the evolution calculations, here restricted to a reasonable 
but uncertain initial value set to Mi = 1.4 Mq for defi- 
niteness. The evolution of the NS (pulsar) and donor star 
masses are depicted in Fig.[3| The final value of the NS mass 
is not sensitive to Pi as expected, and refiects the integrated 
accretion over the system lifetime. While values ^2.2 Mq 
result, it should be kept in mind that the accretion efficiency 
parametrised by /? was held fixed at 0.5 in all the calculation. 
A detailed modelling of the accretion could refine this rough 

•^ This is due to the fact that the larger the Pi the star has a 
higher mean molecular weight fi. The star attains a partially de- 
generate interior at larger densities, a smaller Roche lobe and 
consequently a shorter orbital period. Partial degeneracy is a nec- 
essary conditio n for the orbit of the system to start to ev olve to 
larger periods l lSavoniie. de Kool. fc van den Heuvellll986h . 
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time-average and eventually push the mass to higher values 
without any serious conflict with present observations. How- 
ever, the problem of explaining t he presence of very m assive 
NS in these systems would call (|Romani et al.ll2012|) for a 
major revision of current NS c ore microphysics ( HorvathI 
I2OO6I : iLattimer fc Prakash"2007''). 

In Fig.|4]we show the period evolution for the considered 
CBSs. This figure, complementary of Fig. (T) provides the 
time scale involved in the evolution of the CBSs. Finally, in 
Fig. Owe show the evolution of the radius of the donor star 
and its corresponding Roche lobe for Pi = 0.88 d, which 
represents a good model for PSR J1311-3430. 

In Table [T] we present some relevant quantities for each 
of the computed models. From that, we see that all mod- 
els with Pi < Pb reproduce the main characteristics of 
BW systems. Particularly, the system containing the pul- 
sar PSR J1311-3430 is best represented by models with 
Pi — 0.870 — 0.880 d, just below the bifurcation value. 

Let us briefly address the effect of the present uncertain- 
ties in the ingredients of the model. The value of Pi favoured 
for the occurrence of PSR J1311-3430 is dependent on the 
strength of the magnetic braking. We have tested it by intro- 
ducing a factor lower than one in the standard expression. 
If the strength is half of the standard one, the whole sce- 
nario is essentially unchanged; the favoured period is 0.76 d 
and companions to BWs without H are possible. However 
if magnetic braking were four times less efficient it is still 
possible to find systems with the observed period but in all 
cases the companion retains some H. Thus, if so, a model 
for PSR J1311-3430 would be not possible. Regarding the 
value of the parameter / Lp, work in progress shows that if 
it were larger than assumed here, the qualitative behaviour 
of the evolution would be essentially the same, but the final 
evolution to large orbital periods would be markedly faster 
posing no problem for the whole model. 

Performing a detailed calculation of the probability of 
having BWs without (or trace) H is beyond the scope of 
this paper. However, in order to make some preliminary es- 
timations, we have computed several models for the range 
of periods that evolve to BW configuration. We found that 
for a standard magnetic braking, systems evolving to BW 
have 0.500 ^ P^/d < 0.887; being 0.873 s: Pi/d < 0.887 
for X < lO-'^ (0.865 sj Pi/d < 0.887 for X < 10"^). 
Meanwhile, for the case of magnetic braking with half of 
the standard strength, the range of Pi evolving to BWs 
is 0.500 ^ Pi/d ^ 0.769; being 0.760 ^ Pi/d ^ 0.769 for 
X < lO"'^ (0.753 s; Pi/d ^ 0.769 for X < 10"^). So, if the 
initial distribution of Pi were uniform, we would have 6% of 
the BWs with X < 10"^ and 3% with X < 10"^ for both, 
the standard magnetic braking strength and a half of it. 



3 CONCLUSIONS 

A new path for CBSs evolving into millisecond pulsar- 
very low mass companion final states has been stud- 
ied and applied to the case of P SR J 1719- 1438 in 
iBenvenuto. De Vito. fc HorvathI l|2012r ). In this Letter, we 
have applied the same ideas to the latest member just dis- 
covered in gamma-rays (|Pletsch et al.l 2012'). namely the sys- 
tem PSR J1311-3430. These self-consistent calculations in- 
cluding all three essential ingredients (accretion, winds and 
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Figure 1. Orbital period-mass relation for the 2Mq normal 
donor star evolving together with a 1.4Mq neutron star for dif- 
ferent values of the initial periods (in days). The observed orbital 
periods of PSR J1719-1438 and PSR J1311-3430 are denoted with 
dashed lines. 



illumination) show that theory is able to account for the 
formation of BW systems with He-dominated donor compo- 
sition on a few Gyr timescale. 

The initial conditions for the particular case of 
PSR J1311-3430 are tightly bound from below by the re- 
quirement that X < 2.5 X 10~^, as displayed in Fig. [2l We 
also assert that the same evolution leading the system to the 
observed region of Port — M2 plane produces a high-mass pul- 
sar as indicated by dedicated observations bv lRomani et al.l 
(|2012r ). Figs. IH5I displav the evolution of the system at a 
glance. 

We have finally checked that a variation of the initial 
donor mass M2 within the whole interval leading to the BW 
track does not lead to dramatic changes in the initial pe- 
riods, and that a bifurcation locus can be calculated quite 
accurately. Work in progress devoted to the whole popula- 
tion of "r edbacks", wide binaries and BW millisecond pulsar 
systems ({Robert sll2012l ) is underway and will be published 
elsewhere. 

J.E.H. has been supported by Fapesp (Sao Paulo, 
Brazil) and CNPq, Brazil funding agencies. 
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Table 1. Some relevant characteristics of the binary systems when they reach the observed orbital period for second time. We list the 
initial period, the age, the donor and NS masses, the central temperature and density of the donor star, the orbital period, the surface 
H abundance, the radius of the Roche lobe, the radius and the mean density of the donor star. For the case of initial periods above 
bifurcation we list these quantities at the age of 13 Gyr. 
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Figure 2. The evolution of the H abundance of the donor star 
at its center (thin lines) and surface (thick lines) for the models 
included in Fig. [T] Depending on the value of the initial period, 
at a given age the star becomes fully convective being chemically 
homogeneous. 



Figure 3. The evolution of the masses of the donor star (lower 
panel) and neutron star (upper panel) for the same systems in- 
cluded in Fig. [T] The labels of the different line types denote the 
values of the initial period (in days). 
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Figure 4. Orbital period evolution for the same systems included 
in Fig. [1] 
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Figure 5. The radius of the donor star R2 (solid lines) and its 
corresponding Roche lobe Rj^ (dashed lines) for the case of an 
initial period of 0.88 days. Points denote the conditions at which 
the system attains the orbital period of PSR J1311-3430 with low 
mass in detached conditions. 
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